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1. Introduction 

Recent technological advances make it possible to realise simple mechanical devices in 
the microscopic and nanoscopic regimes, the properties of which are determined by 
quantum effects [T]. The existence of these represents a remarkable opportunity for 
fundamental studies of light-matter interactions [2] and also the potential for practical 
application to quantum communications and information processing [3]. Yet they 
present also a challenge to existing methods of analysis, many of which were developed 
to treat more rapidly oscillating systems with weaker couplings. 

The quantum theory of machines is built, to a large extent on the theory 
of oscillators and strongly coupled oscillators, which are coupled to one or more 
environments each of which is at a characteristic temperature in El E]- The behaviour 
of these quantum systems is governed not simply by average properties but also 
fluctuations, and has been informed by the development of fluctuation theorems for 
quantum open systems HIE]. The coupling to the environment needs to be handled 
with some care, however, because of the possibility of quantum coherence and the 
development of entanglement between the system of interest and its environment, with 
the result that apparently unphysical behaviours may emerge [9|. The requirement to 
treat the coupling to the reservoirs with care provides the incentive to return to the 
problem of a single strongly-damped oscillator and to treat this model exactly. 

There exist at present mathematically and physically acceptable methods for 
treating strongly damped quantum oscillators, but most of these do not retain the 
simplicity of their classical counterparts. Our aim is to recover at least some of this 
simplicity. We do not seek to apply familiar techniques developed for weakly damped 
systems, but proceed instead by means of an exact diagonlisation of the oscillator- 
reservoir Hamiltonian. This approach is complementary to that adopted by Philbin who 
has tackled this problem of the oscillator evolution using Green functions [Tn|. We hope 
that the combination of his work and ours will provide a more complete understanding. 

Our procedure produces a single positive function of frequency the integral of which 
is unity. We can interpret this quantity as the frequency spectrum or probability density 
for the oscillator in its ground state. The steady-state of the oscillator, its entanglement 
with the environment, its energy and much of its dynamics can be understood in terms 
of this probability. Before embarking on this programme, let us pause to review some 
of the background to our study. 

1.1. Background 

The harmonic oscillator has a special place in physics as one of the simplest and most 
widely employed of physical models. The reasons for its ubiquity, no doubt, are its 
simplicity and the fact that it is readily analysed. In the quantum domain, the harmonic 
oscillator is barely more difficult to treat than its classical counterpart and was one of 
the hrst dynamical systems to which Schrodinger applied his equation m- Today, both 
the classical and quantum forms appear in elementary courses on classical and quantum 


Theory of the strongly-damped quantum harmonic oscillator 


3 


mechanics. 

The damped harmonic oscillator loses energy as a result of coupling to the 
surrounding environment. In the classical domain it often suffices to describe this in 
terms of a simple damping coefficient, 7 , and an associated stochastic or Langevin 
force [ 12 ], ^{t), which models the effect of environmental fluctuations on the oscillator 
[la El na uBi HD. The dynamics is described by a simple linear differential equation 
of the form 

x + 'jx-\-uix= -, (1) 

m 

where m is the mass of the oscillating particle. There is no requirement for detailed 
knowledge of the fluctuating force, which may be considered to have a very short 
correlation time with a magnitude determined by the requirements of thermodynamic 
equilibrium. 

The damped quantum harmonic oscillator requires that explicit account be taken 
of the quantum nature of the environmental degrees of freedom [TH], which are most 
simply described by an ensemble of harmonic oscillators [19]. If the damping is very 
weak, so that 7 wq, then we can neglect rapidly oscillating terms in the coupling 
between the oscillator and the environmental oscillators by making the rotating wave 
approximation, which corresponds to enforcing the conservation of the total number of 
vibrational quanta, and then the Born and Markov approximations associated with weak 
coupling and loss of memory in the reservoir [20]. This leads to the master equations 
and Heisenberg-Langevin equations that are ubiquitous, most especially, in quantum 
optics [iniEIlEllEalElESlEniETlEElEHllaoll^ 

If the coupling is somewhat stronger then it may not be possible to make the 
rotating wave approximation and we then need to retain in the Hamiltonian terms 
that can create or annihilate a pair of quanta, one in the damped oscillator and, at 
the same time, one in the environment. This leads to the Caldeira-Leggett model 
[In] [32] l33] IMl l35] l36] [37], which we describe in the the following section, and which 
has been applied to study a wide variety of quantum open systems [271 El]. A variant 
on the model has been applied to the quantum theory of light in dielectric [3HI [39] and 
magneto-dielectric media An important complication that seems to 

be an inevitable consequence working in this strong-coupling regime is the failure of 
the Markov approximation; attempts to enforce this approximation lead to a master 
equation that is unphysical in that there exist initial states for which the dynamics 
leads to negative probabilities [271 01 ig 061 ST]. It is possible to derive a master 
equation but the resulting equation is one that has within it non-trivial time-dependent 
coefficients [371 sg. This time-dependence is a clear signature of the non-Markovian 
nature of the associated evolution. It seems that this non-Markovian character is an 
inevitable feature of the strongly-damped quantum harmonic oscillator. 
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Figure 1. Representation of a single harmonic oscillator (in red) coupled harmonically 
to a bath of oscillators (in blue). 


2. Hamiltonian for the strongly-damped harmonic oscillator 


Consider a harmonic oscillator of natural angular frequency uo that is strongly coupled 
to a large collection of oscillators, with a range of frequencies, through their respective 
positions as depicted in hgure [T] We write the Hamiltonian for the combined oscillator- 
reservoir system in the form 


H = ——h -mu^x^ -b 
2 m 2 ° 


^ Pi 

2m, 




- Y, m^ulXf,Xf,x. 


( 2 ) 


If we complete the square we can rewrite this in a minimal-coupling form to arrive at 
the alternative form: 




Pi 

2m„ 


+ 


{x^ - X^xY 


where 


o2 - , ,2 


E Aj- 


m 


( 3 ) 


( 4 ) 


Each term in our Hamiltonian is strictly positive only if this quantity is positive. If it 
is negative then the second term is also negative and the Hamiltonian is not bounded 
from below and hence not allowed physically. Hence the positivity of the Hamiltonian 
places a physical restriction on the strength of the coupling 
2 ^ 2 \2 


m 


( 5 ) 






Theory of the strongly-damped quantum harmonic oscillator 


5 


It is convenient to rewrite our Hamiltonian in terms of the familiar annihilation 
and creation operators: 

“ ^ ip \ 


a = 


= 


moof) 


X + 


mojQ 


2h 


Xn + 




( 6 ) 


In terms of these operators our Hamiltonian becomes 

H = hujod'^d + (« + {K + K) 

ti fi ^ 

when unimportant constant shifts in the ground-state energies are removed and 


n = 


mooo 


( 8 ) 


( 9 ) 


When written in terms of this quantity, our positivity condition (|^ becomes 

1/2 

At ^At 

There is one further rehnement that we need to introduce before considering the 
diagonalisation of the Hamiltonian and this it to introduce a continuum description 
of the environment. To this end we introduce continuum annihilation and creation 
operators, b{u;) and b^{u;), satisfying the commutation relations 

b{uj),b\u') = S{uj — u') (10) 

and our Hamiltonian becomes [39] 


( 11 ) 


/V /*CO 

H = huoofd -f / doj hoj¥{ijj)b{pj) 

Jo 

POO , \ ^ 

J du — V (cn) (a -1- a^j \b\uj) -|- b{uj) 


and the positivity condition is 


uJo> duj- 

Jo OJ 


( 12 ) 


Our Hamiltonian is quadratic in the annihilation and creation operators for the oscillator 
and the reservoir and hence leads to linear coupled equations of motion for these 
operators. We could seek to solve these equations of motion and this would lead to an 
operator Heisenberg-Langevin equation analogous to ([^. Here we adopt the different 
approach of diagonalising the Hamiltonian. 


2.1. Hamiltonian diagonalisation 

We seek to diagonalise the Hamiltonian by hnding a complete set of eigenoperators, 
B{uj) and their conjugates i?'^(a;) that satisfy the operator equations 

B (u ), H = huB (u) 

'b\uj),h] =-huB^iuj) (13) 
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for all positive frequencies u. These eigenoperators are complete if, in addition to these 
they also satisfy the condition 

'b{u),B\u')]=6{u-u'). (14) 

These operator equations are the natural analogues of the more familiar eigenvalue 
and completeness conditions for the eigenstates of a Hamiltonian. In analogy with 
the eigenvalue problem, we expand each of the eigenoperators as a superposition of a 
complete set of operators: 

B{u)) = a{u)d-\-(3{u)af / dto' 5{u^u')b\u') , (15) 

Jo L 

and then use the eigenoperator equations and completeness condition to determine the 
coefficients in this expansion. The calculation is a little involved but the main points 
are summarised in [Appendix A 


We can express any of the annihilation and creation operators for the oscillator or 
the reservoir in terms of our eigenoperators. To do this we write the desired operator 
as a superposition of all the B{u) and i?l(a;) operators and then use the commutation 
relations to extract the coefficients in this expansion. For the oscillator operators we 
hnd 


a = 


d^ = 


dto (a*{u)B{u) — (3 {u)B''{u 
du [a{u)B^{u) — (3*{u)B{u)j . 


(16) 


The requirement that these operators satisfy the familiar boson commutation relation. 


a, a 


= 1, provides a constraint on the functions a{u) and /3{uj) in the form 

dUoCJ 


du 


\a{u)f - \P{u)Y 


du\a{u)Y 


= 1 , 


(17) 


(cuo + uY 

where we have used ( A.10| ). It is interesting to note that the correctness of this may be 
verihed explicitly by contour integration, but that the proof makes explicit use of the 
positivity condition ( [I^ [39]. 

The integrand in ( [l7j ) is clearly positive (or zero) for all frequencies and hence has 
the mathematical form of a frequency probability distribution: 

AuqU 


7r(a;) = |a(c(;)|" 


( 18 ) 


A number of further constraints on this quantity emerge naturally from thinking of it 
as a probability distribution and from our diagonalisation. 


2.2. Physical constraints 

We leave until the next section the question of the physical interpretation of the 
probability density vr(a;) but consider here what can be inferred from the fact that 
it has the mathematical properties of a probability distribution. To this end let us 
denote the average value of a function of frequency for this distribution by 


roo 

((/(^))) = / duf{u)7i{u). 
Jo 


( 19 ) 
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we note that |Q!(a;)p is finite for all u and it follows, therefore, from (18) that the average 
value {{u~^)) is finite. 

It follows from the eigenoperator equations (13) that we can write the Hamiltonian 
in the form 

.V poo 


H = 


10 


du hu w{uj)B{uj) + C, 


( 20 ) 


where C is an unimportant constant. If we substitute into this Hamiltonian our 
expression for the eigenoperators (15) we require that the coefficients of of and 
should vanish so that 


/ du u\a(uj)\^ 
Jo 

which implies that 

{{u^))=ul 


2 fU — Uq 
U + Uq 


, niu), 9 9 , 

du— - iu^ — Un) = 0 

dwn 


( 21 ) 


( 22 ) 


The fact that the square of the mean value cannot exceed the mean value of the square 
for any probability distribution leads us to deduce that 

((a;))<a;o. (23) 

Finally we can apply the Cauchy-Schwartz inequality to provide a lower bound on 
the value of {{u~^)): 

{{u-^)){{u)) > 1 


Uq 

These inequalities are useful in determining the properties of the ground state. 


(24) 


3. Ground-state 

It is immediately clear from the form of the Hamiltonian ( [IT| ) that the ground state of 
the oscillator is not that of the undamped oscillator, that is the state annihilated by a. 
To see this we need only note that the interaction term a^b^u) acts on the combined 
ground state of the non-interacting oscillator and reservoir, adding a quantum to each. 
This suggests that the true ground state should be a superposition of states with varying 
numbers of quanta in both the oscillator and the reservoir and hence an entangled state. 
This situation is reminiscent of the ground state of an atom in quantum electrodynamics, 
which is dressed by virtual photons |19]. The dressing of the ground-state atom is 
responsible for a number of important effects including the Casimir-Polder interaction 
[501 El] and the form of the polarizability of the atom [52l E3] . It is reasonable to expect 
that it will be similarly significant for our strongly damped oscillator. 

The true ground state, which we denote by the ket |0), is the zero-eigenvalue 
eigenstate of all the annihilation operators B{u): 

.B(a;)| 0 ) = 0 \/ u. 


(25) 
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The ground-state properties of the oscillator in this pure state are described by a mixed 
state density operator obtained by tracing out the environmental degrees of freedom: 

Posc = TrEnv(|0)(0|). (26) 

The most straightforward way to determine the form of this mixed state is to use the 
characteristic function 


(27) 


X(0 = Tr [pexp 

= {0|exp (^a^' -Ca) |0). 


This function provides a complete description of the state and all of its statistical 
properties. To evaluate it we express a and in terms of the eigenoperator^ 


X(0 = {0|exp 
= exp 


du 


{(a + CP) BHuj) - (Ca’ + ?/ 3 *) \ | 0 ) 




(28) 


This simple form, Gaussian in is characteristic of a squeezed thermal state. We can 
rewrite this characteristic function in terms of our probability density in the form 


(29) 






w 

exp 

-- / 7r(u;) 

—c 

+ — c 



L 2 Jo 

Vwo 

OJ 

J \ 


= exp 


1 

2 




Uo 


-er + {{co-^))u;o^t 


where and are the real and imaginary parts of ^ respectively. We note, in passing, 
that the quadrature operators for the oscillator (familiar from quantum optics [20]) have 
unequal uncertainties: 


A 


A 


a -|- a’*'\ 

' —i{d — al) 

; V 2 ", 


'((a;-i))a;o 


'M), 

2 ujq 


(30) 


The product of these two exceeds as it must, by virtue of the Cauchy-Schwartz 
inequality (24). 

We have established that in the ground state, the reduced state of the oscillator is 
Gaussian in position and momentum and, as such, is completely determined by the first 
and second moments of the position and momentum. These moments are 

{ 0 |x| 0 ) = 0 

(0|p|0) = 0 


(Olf^lO) = 




2m 


f We make use of the operator ordering theorem pO] 

exp(T B) = exp(T) exp(.B) exp B] 

which holds if the two operators A and B both commute with their commutator \A, B\. 
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(OIp^IO) 


hm{{oj)) 

2 


{0|pa; + px\0) = 0. 


(31) 


We note that the frequency of the undamped oscillator, uq, makes no explicit appearance 
in these expectation values. The mean energy of the oscillator alone is 

(0|^ + = ^huo (32) 

The fact that this exceeds the ground-state energy of the undamped oscillator is a 
reflection of the fact that there is an energy cost to be paid in order to decouple the 
oscillator from its environment [9]. We note, further, that the mean kinetic energy and 
potential energy for the oscillator alone do not have the same value and that this is in 
marked contrast to the ground state of the undamped oscillator. 

We can diagonalise the density operator for the oscillator alone, posc, by means of 
a squeezing transformation [20] or, equivalently, introducing a new pair of annihilation 
and creation operators corresponding to a new oscillation frequency, Uc'- 


c = 


= 


moOc I 

~W \ 

muc 


X i- 


P 


mUr 


X 


. p 

(- 

mujc, 


(33) 


To complete the diagonalisation we need only choose ujc such that the expectation values 
of and are zero: 

^2 


moOc 


10 ) = 


2h 


\x 


P 


| 0 ) = 0 


bJr = 




(34) 


This frequency is, by virtue of (23) and (24), less than cuo, the natural frequency of the 
undamped oscillator. The appearance of a reduced oscillation frequency is reminiscent 
of the observed reduction in the oscillation frequency of a damped classical oscillator but 
unlike in the classical case, the frequency uJc remains real even in the heavily damped 
regime. We note that Philbin and Horsley arrived at a similar reduced but always real 
frequency, which they associated with a decrease in the zero point energy 
The mean number of c-quanta in the oscillator ground-state is 

1 


n, 


= {o|ct«|o) = i [fTmTT) 


(36) 


which exceeds 0, as it should, by virtue of (24). When written in terms of the c-quanta, 
the steady-state density operator takes the form of a thermal Bose-Einstein state, which 
we can write in the form 


POsc — _ 


Ur 


Uc+l XUn + 1 


(36) 
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We may interpret this state as a thermal state for the oscillator at the shifted (and 
reduced) frequency Uc and at an effective “temperature” 

We should note, however, that the true temperature in the ground-state is zero and that 
this quantity and the frequency ujc are at most only parameters with which to quantify 
the state of the oscillator and its entanglement with the environment. In particular, the 
von-Neumann entropy associated with the steady state of the oscillator is 

S'(Osc) = (he -f 1) ln(hc -|- 1) — he In he, (38) 

and that, by virtue of the Araki-Lieb inequality [SSI EHl EZj and the fact that the full 
state is pure, means that this is also the total entropy of the environment: 

^(Env) = ^(Osc) (39) 

and that the quantum mutual information ra. or index of correlation jSS], between the 
oscillator and the environment is 

5(Osc : Env) = 5(Env) + ^(Osc) - ^(Osc, Env) = 25(Osc). (40) 


4. Physical meaning of vr(ci;) 

We have seen that the physical properties of the oscillator ground state may readily be 
expressed in terms of moments of the frequency given the probability distribution 7r{uj). 
Here we present the case for identifying this probability density with the contribution 
from the dressed modes, associated with the eigenoperators B{u), to the the state of the 
oscillator. It is in, in essence, the spectrum of the true ground-state modes contributing 
to the oscillator state. We present three arguments to support this interpretation. 


Our Erst justiheation arises from the form of the expectation values (31). We know 
that the ground state of a harmonic oscillator of frequency u has 

h 

2mu 
hmu 


{x^) = 

{f) = 


(41) 


If we treat the state of the oscillator as a mixture of oscillators of different frequencies, 
each in its ground state, and contribution with weight vr(a;) then the resulting average 
mean-square values will be 

h h{{u~^)) 


(x^) = / duj7r(uj) 

^ ^ Jo ^ ’2muj 

,.2\ /■“ j / .hmu 

{p)= du vr(a;) 

Jo 


2m 

hm{{u)) 


10 2 2 ’ 
which correspond to those obtained above. 

Our second point arises from the form of the Hamiltonian for the oscillator 


( 42 ) 
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We can, by virtue of (22), write this as a combination of potentials corresponding to 


different frequencies but weighted by 7r{u): 


H„=^+hn 


du 7r{uj)u}‘^x‘^ 


2m 

p2 

2m 


-m{{u'^))x\ 


(44) 


Finally we note that the effective mean energy of the oscillator, which is associated 
with the diagonal form of the density operator ( p^ is 

1 \ . 1 . /■“ . . . 1 , 


n, 


1 \ 1 1 


(45) 


which combines the characteristic ground state energies of the dressed oscillators, 
weighted by the probability distribution 7r(a;). The combination of these three features 
leads us to interpret vr(a;) as the proportion of the corresponding dressed oscillators 
contributing to the properties of the damped oscillator. We emphasise that the 
mathematical results obtained in the preceding section do not require us to adopt this 
interpretation of vr(c(;) but we find it helpful to do so. 


5. Oscillator dynamics 

The exact diagonalisation of the Hamiltonian makes it straightforward to evaluate the 
time-evolution of any desired property of the oscillator. All that we need do is to express 
the desired observable in terms of the eigenoperators, B{cj) and and then use the 

time evolution of these operators, the form of which is an elementary consequence of 
the fact that they are eigenoperators: 

B{uj,t) = 

W{u,t) = B^{u,0)e^‘^h (46) 

In particular, we can determine the time-evolution of the annihilation operator for the 
oscillator in this way: 


iujt 


d(t) = / du a*{u)B{u,0)e — /3{u)B^{u,0)e''' 

Jo L 

/*OCi f p 

= J dcaj a*(a;)e“*^* a(c(;)a(0)-|-/d(a;)a^(0)-|- 

J du' (pf{u,u')b{u',0)6{u,u')b\u',0)'^ 
— /3(a;)e*‘^* a*(a;)a^(0) -|- /?*(a;)a(0)-|- 

poo / ^ 

/ du' (pf*{u,u')b\u',0)6*{u,u')b{u',0) 


(47) 


This may be used, together with the initial state of the oscillator and the environment, to 
evaluate the expectation value of any desired property. We note that if the environment 
is in a stationary state, so that (6(a;,0)) = 0 = (6l(a;,0)), then the expectation values 









Theory of the strongly-damped quantum harmonic oscillator 


12 


of the position and momentnm operators take a pleasingly simple form: 

{x{t)) = ((cos(a;f)))(f(0)) + — ((o;"^ sin(a;f)))(p(0)) 

m 

{p{t)) = {{cos{ut))){p{0)) - m{{usm{ut))){x{0)), (48) 


where the donble angle brackets denote, as before, averages over onr probability 
distribntion vr(a;) as in (19). The form of these eqnations adds farther snpport to onr 
interpretation of vr(a;) as a freqnency probability distribntion for the damped oscillator, 
as they may be viewed as the evolution of an undamped oscillator with a frequency 
uj averaged using this probability distribution. The dissipation arises simply from a 
dephasing amongst the different frequency components. 

The evolution of the mean position and momentum, as given in (48), has the 
necessary short-time form of that for an undamped oscillator 


{x{dt)) = (h(0)) + —(p(0)) 
m 


(49) 


{p(dt)) = {p(0)) - muildt{x(a)), 

where we have used the identity ((ca^)) = Uq. The effects of the coupling to the 
environment enter at order dt^ and this is an indication of the essentially non-Markovian 
nature of the strongly-damped oscillator. Our interest is in strongly damped oscillators 
and so we should note that (48) includes the possibilities of both critically-damped 
and over-damped evolution. The equations contain, moreover, a simple criterion for 
these, which we may express in terms of our probability distribution. The motion 
will be oscillatory if {{cos{ut))) has stationary points at times other than at t = 0. 
Alternatively, we may state that the motion is under-damped if the derivative of this 
quantity, that is ((ci;sin(ci;f))), is zero for any time other than t = 0. If it is zero only at 
t = 0 then the motion is critically-damped or over-damped. 

Our expression for the evolved annihilation operator (47), together with the 
corresponding one for the creation operator provide a full description of the oscillator 
dynamics. This is true for any initial state of the oscillator and, moreover, for any 
environmental state including, of course, that associated with a hnite temperature. As 
an illustration let us examine the evolution of the characteristic function for an arbitrary 
initial state of the oscillator coupled to a zero-temperature environment at time t = 0. 
With a little effort we hud 


1 

roo 

2' 

-- / du' 

/ dco fi(u, u', t) 


2 Jo 

Jo 



where 


X(CT) = X[^(^),0] exp 


^{t)= / du[ri{uj,t)a*{uj)-ri*{uj,t)f3{uj)] 

Jo 

oj\t) = r]{u, t)'y* (cj, u') — p* {u, t)6* {u, u') 

7^(a;,t) = [e«(a;)+r/5(a;)]e*‘^h 


(50) 


(51) 


This characteristic function encodes the full dynamics and statistics of the oscillator. 
As a simple illustration of this we can determine the form of the steady state. To see 
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this we first note ^{t) tends to zero as t tends to infinity and the different freqnency 
components dephase so that 

X[e(oo),0]=x[0,0] = l. (52) 


Evalnating the long-time limit of the exponential factor in (50) reqnires some care in 

)art components. We £n 

den |^a(a;') + (53) 


the handling of the delta-fnnction and principal part components. We find 


dco' 


du bj\ oo) 


so the steady-state characteristic fnnetion is 

1 


xK.oo) =exp (“ 2 ^ dhj\(a(hj) +r/3(w)|^) , 


(54) 


which we recognise as the characteristic fnnetion for the oscillator in the global ground 


state (28). This is a most satisfactory and exact result. 


6. Weak-coupling limit 

The theory developed above was designed to treat the strongly-damped harmonic 
oscillator, but should also be applicable to the more familiar weakly damped oscillator, 
for which the oft-employed Born and Markov approximations are applicable and the 
steady state of the oscillator should be its ground state. We show here that this is 
indeed the case. 

We start by considering the form of the function a{uj) in the weak coupling limit. 


To aid our analysis we rewrite the form given in (A.18) as 

U!-\-U!o f V{U!) \ 

a{oj) = ' ' 


(55) 


uq \\V {u)\'^Y (u) — iTr\V J 

The weak damping limit corresponds to choosing the coupling to the environment to be 
small or, more specifically, to |l/(a;)p <C ujq. It is clear that in this limit, |q;(ci;)P will 
be a sharply peaked function centred around the frequency for which Y (cu) = 0. If the 


integral part in Y(u), as given in (A.15) is smal|§]then this frequency will be close to 
Uo and we can write 

\V{uj)fY{uj) ^ 4{u - ooo) - 4F{uj), 


where 


4F{u) = den' - 
Jo \oj 


p 


■ U' 




(56) 


(67) 


(58) 


UJ -\- Ujf 

This leads, in turn, to a corresponding approximate form for a{uj): 

. 5 V{u) 1 

2 a;-cno-F(cn)-^f|E(cn)r 
where we have set u = ojq everywhere except in the denominator. We note that this 
is of the form that arises from the Fano diagonalisation of our problem if we make the 

§ If it is not then we will need to invoke ideas of renormalisation, a complication we wish to avoid. 
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rotating wave approximation by omitting from our original Hamiltonian all terms that 
are products of two creation operators or of two annihilation operarors [20l l59] . 

Consistency with the above approximation, which led us to set oj = ojq leads us to 
set (3{oj) to zero: 


P{u) ^ 0 , 


(59) 


so that the integral over all frequency of \a{u)\‘^ is unity. Hence in this limit the steady- 
state characteristic function for our oscillator is 



which we recognise as that for the ground state of the undamped oscillator [ 20 ] . 
Further, we note that in this limit our probability distribution function, 7 r(a;) ~ \a{u)\‘^, 
approaches a Lorentzian centred on uq, with some width 7 [j|| Thus all the complexity 
of of the original problem is reduced, in the weak-damping limit to just two parameters: 
a natural oscillation frequency and a damping rate. 

It was important to conhrm that our more general treatment coincided, in the right 
limit, with the approximate methods used for weakly damped oscillators. We should 
note that even if we are working in the weakly damped regime, then our approach offers 
a systematic way to treat corrections to the results obtained using the Born and Markov 
approximations, which may play an important role in modelling measurements at the 
limits of sensitivity. 

7. Conclusion 

We have presented an exact diagonalisation of a simple quantum model of the damped 
harmonic oscillator, one that is applicable, in particular, to strong damping. We have 
found that much of the behaviour of the oscillator and may of its properties can be 
described in terms of a single probability function, 7r(a;), which we may interpret 
as the contribution of corresponding dressed mode, at frequency u to the oscillator. 
These properties include the steady state of the oscillator at zero temperature, the 
entanglement between the oscillator and its environment and also its evolution, both 
in the familiar under-damped regime but also in the more problematic over-damped 
regime. 

We have applied our diagonalisation to study the properties of the true ground state 
and have shown that the oscillator part of this pure entangled state coincides with the 
steady-state of the oscillator in a zero-temperature environment. The diagonalisation 
is not specihc to any particular state of the reservoir, however, and may readily be 
applied to environments at hnite temperature. It may be extended, moreover, to include 
driving forces, coupled oscillators and multiple reservoirs, with the latter perhaps being 
at different temperatures [60]. This may provide some insights into important questions 


II This is not strictly true in the wings of the distribution, of course, as even in this limit we require 
((w^)) = ujq, but the corresponding quantity for a true Lorentzian is divergent. 
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of principle in the nascent fields of quantum machines and quantum thermodynamics 
[611 [621163| IMl l65l [66] . We return to these ideas elsewhere. 
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Appendix A. Fano diagonalisation 


We present a brief account of the exact diagonalisation of our Hamiltonian based on 
methods developed by Fano for the study of conhguration interactions [HZj. This idea 
was applied extended to weakly-coupled oscillators in a quantum study of damped cavity 
modes [20l [59]. The extension to stronger couplings, with the inclusion of counter¬ 
rotating couplings, has been given before and applied to the quantum theory of light in 
dielectric media [38l [39] • We summarise here the analysis presented in [39| . 

Our task is to diagonalise the damped harmonic oscillator Hamiltoniaij^ 


/V /*00 

H = huoa^d du hub\u)b{uj) 
Jo 


+ 


du — (jiV{u)P{u) + V*{u)b{u) 


(A.l) 


by which we mean rewriting it in the form of a continuum of uncoupled or dressed 
oscillators: 


/V /*00 

H= / duhu BUu)B(u)+ C, (A.2) 

Jo 

where C is an unimportant constant. We proceed by writing the dressed annihilation 
operators, B{u), as linear combinations of the bare operators for the oscillator and bath 
modes: 


B{u) = a{u)dft{u)d ^/ du' 'y{u,u')b{u') 6{u,u')¥{u') , (A.3) 


JO ■- -■ 

where au, (J{u), 'y(u,u') and 6(u,u') are to be determined. 

We require the operator B(u) to be associated with an uncoupled or dressed 
oscillator of angular frequency u. This requires us to hnd its form such that the following 
pair of operator equations are satished for every frequency, u: 


B {u ), H = huB (u) 
B{u),B'^{u') =6{u — u'). 


(A.4) 

(A.5) 


^ For the Hamiltonian of interest in this paper the coupling, F(a;), is real but treating the problem 
with a more general complex coupling presents no additional difficulties. 










Theory of the strongly-damped quantum harmonic oscillator 


16 


Substituting the ansatz (A.3) into (A.4) and comparing coefficients of the bare creation 
and annihilation operators leads to the set of couple equations: 

X 


1 

a{ijj)ojQ + - du' = a{u:)ijj (A. 6 ) 

2 JO 

—I3{u:)u:q H— / du' [^{u,u')V{u') — 5{u,u')V*{u')] = ld{u)u (A.7) 
2 Jo 

a{u) — ld{u)] + 7 ( 0 ;, u')u' = 7 (a;,a;')a; (A. 8 ) 

a{u) — jdiy)] — 6(u,u')u' = 6{u,u')u. (A. 9 ) 


V*{u') 

2 

V{u') 


Our method of solution is use these to determine the functions fdu, 7 ( 0 ;, cj') and 
S{u,u') in terms of au and then to determine this remaining function by enforcing 
the commutation relation (A.5). From (A. 6 ) and (A.7) we see that 


/3{u) = -^ a{u). 

^ ^ CJ + CJo 

If we use this to substitute for (3{u) into the remaining equations then we hnd 


V*{u') 

V{u') 


U + Uq 
Uq 


a{u) = 'y{u,u'){u — u') 


(A.IO) 

(A.ll) 
(A. 12 ) 

. . (A.13) 

CJ + CJo 

The hrst, however, requires careful handling because the behaviour at cj = cj'. Following 
Fano isa. we adopt the method proposed by Dirac [SB] and write 


-Q;(a;) = 5(a;, u'){u + u'). 


CJ + CJo 

Solving the second of these presents no difficulty and we hnd 

5{u,u') = ( -^ - - — a{u). 

^ ’ \u + u') ^ ’u + uq ^ ’ 


7 ( 0 ;,cj') = [-- + F(a;)(5(a; — cj') J — 0 !{u), (A.14) 

where P denotes that the principal part is to be taken on integration and Y [u) is a real 
function, which we determine by substituting (A.14) into (A. 6 ). We hnd 

1 


Y{u) = 




2(a;2-u;2) r 


du' 


P 


Uo 


u — u' U -\- U‘ 


\VU)\ 


(A.15) 


If we substitute our operators, B{u), expressed in terms of the function Q;(a;) into 


the commutation relation (A.6) then we hnd 

B{u),B^{u') = a;(a;)Q:*(a;') |l— 
P 


CJ — CJo \ / Cj — CJo 


+ 


du 


^u — U‘ 
1 


+ Y{u)6{u- u") 


U -\- UqJ \u' + Uq 
p 


u' - u" 


+ Y(u)S(u' - u") 


u + u" J Vcj' + U' 


\V{u")\^ul 


(cu + a;o)(a;' + a;o) j ' 


(A.16) 
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Evaluating the integrals and setting the result equal to (5(a; — ca') give^ 

I f na _ \ 


(A.17) 


Note that the diagonalisation does not £x the phase of the complex function a{u) and 
we are free to choose this as we wish. A convenient choice is to set 

1 \ 


Q!(a;) = 


ui 


UlQ 


IjJqV*{oj) \Y{oj)—i'K^ 


(A. 18) 
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